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In  this  project,  we  explored  numerical  methods  for  solving  multi-dimensional  deconvolution  and 
inverse  scattering  problems,  and  were  particularly  interested  in  techniques  based  on  the  singular 
value  decomposition  (SVD);  our  emphasis  was  on  both  speed  and  accuracy.  We  considered  new 
techniques  for  numerical  multilinear  algebra,  in  particular  for  computing  the  SVD  of  a  compound 
matrix,  and  investigated  implementation  on  parallel  hardware. 

We  developed  a  new  algorithm  for  solving  the  boundary  distortion  problem  in  image  processing. 
Image  restoration  is  an  ill-conditioned  process  concerned  with  recovering  an  unknown  scene  from 
an  observed  image.  Much  effort  has  been  devoted  to  finding  robust  (against  noise)  methods  that 
maintain  certain  desirable  properties  such  as  sharpness  in  the  restored  image.  However,  efficient 
restoration  techniques  often  introduce  errors  in  the  modeling  of  image  boundaries  and  these  errors 
are  then  grossly  magnified.  There  are  few  known  methods  for  reducing  boundary  distortion.  We 
invented  a  new  technique  called  incomplete  deblurring  that  modifies  the  assumed  blur  near  the 
boundaries.  Our  approach  is  novel  in  that  we  use  local  (instead  of  global)  regularization,  thus 
getting  the  desired  result  without  excessively  blurring  the  image  in  the  middle  region.  In  [2]  and 
[3],  we  presented  experimental  results  comparing  our  technique  against  other  methods,  and  used 
the  singular  value  decomposition  to  explain  why  our  technique  gives  superior  results. 

This  project  supported  one  PhD  student,  David  Vandevoorde. 
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